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tron  beam  was  used  to  probe  the  plasma  created  by  a  laser  beam  impinging  on  a 
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SUMMARY 


This  report  contains  the  results  of  experimental  work  performed  under  Contract 
DA-19-020-AMI-00521(X) ,  "Laser  Induced  Surface  Effects."  A  moderate  energy 
electron  beam  was  used  to  probe  the  plasma  created  by  a  laser  beam  Impinging  on 
a  graphite  target.  The  attenuation  of  the  beam  was  used  to  determine  the  amount 
of  material  traversed  by  the  beam.  Although  the  original  expectation  was  that 
this  technique  would  be  extremely  useful  in  probing  the  vapor  near  the  target  at 
times  which  were  comparable  to  the  overall  laser  pulse  duration,  the  results  in¬ 
dicate  that  the  bulk  of  the  observed  results  were  due  to  vapor  thermally  emitted 
from  the.  target  at  times  which  were  long  compared  to  the  laser  pulse.  This  effect 
has  been  observed  for  graphite.  However,  no  experiments  have  been  performed  which 
could  confirm  this  effect  for  other  materials.  The  experimental  work  on  the  elec¬ 
tron  beam  technique  was  initially  performed  by  Dr.  H.  Furumotu.  The  experimental 
results  reported  herein  were  obtained  by  L.  Pettingill  and  J.  Shumsky.  The  inter¬ 
pretation  of  the  results  and  preparation  of  the  report  was  performed  by  R.  Schlier. 


I. 


INTRODUCTION 


During  the  past  several  years,  a  number  of  different  techniques  have  been  used  at 
Avco  to  Investigate  the  physical  phenomena  which  occur  during  and  after  the  inter¬ 
action  of  intense  laser  radiation  with  solid  targets.  All  of  these  techniques 
have  been  intended  to  probe  the  properties  of  the  vapor  created  by  the  interaction, 
since  the  vapor  dominates  the  interaction  at  high  intensity  levels. 

The  experiments  that  have  been  performed  include:  (1)  measurement  of  impulse 
transferred  to  the  target;  (2)  measurement  of  velocity  distribution  of  the  plasma 
created  by  the  interaction,  at  times  long  when  compared  to  the  laser  pulse;  (3) 
measurement  of  electron  density  in  the  plasma  during  and  after  the  laser  pulse; 
and  (A)  measurement  of  electron  temperature  in  the  plasma.  In  the  present  report, 
a  fifth  measurement  technique  is  described  —  this  consists  of  measurement  of  mass 
density  of  the  vapor  at  distances  close  to  the  target. 

Before  describing  the  actual  experiment,  it  would  be  well  to  state  the  understanding 
of  the  phenomena  as  it  existed  at  the  time  the  experiment  was  undertaken.  Let  us 
consider  a  solid  target  which  is  irradiated  by  a  short-duration,  intense  laser 
beam.  If  the  target  is  metallic,  the  laser  radiation  is  reflected  and  absorbed. 

Since  the  metallic  target  has  a  large  extinction  coefficient  for  optical  radiation, 
it  is  reasonable  to  assume  that  the  absorption  of  the  radiation  occurs  at  the 
geometrical  surface,  rather  than  in  the  volume  near  the  surface;  consequently, 
the  surface  is  rapidly  heated  until  appreciable  evaporation  can  occur.  However, 
the  temperatures  that  can  be  attained  by  the  surface  are  sufficiently  high  that  a 
large  amount  of  thermionic  electrons  and  metallic  ions  are  emitted  from  the  hot 
surface.  These  ions  and  electrons  fora  the  basis  for  further  absorption  of  laser 
energy  by  the  vapor.  In  particular,  the  electrons  present  can  gain  sufficient 
energy  from  direct  excitation  by  the  intense  radiation,  causing  further  ionisation 
to  occur.  Tne  vapor  then  becomes  highly  ionized  and  very  absorbing  to  the  radia¬ 
tion.  Thus,  the  larger  part  of  the  laser  energy  is  deposited  in  the  vapor,  rather 
than  the  solid,  and  the  vapor  attains  a  very  high  tamperature  and  a  very  high  ex¬ 
pansion  velocity. 

This  qualitative  description  has  been  based  on  the  experimental  observations  ^ 
mentioned  above.  The  expansion  velocities  observed  are  on  the  order  of  3  x  10  cm/sec, 
corresponding  to  particle  energies  of  several  hundred  electron  volts.  The  observed 
vapor  temperatures  are  as  high  as  AO  electron  volts.  The  electron  density  mea¬ 
surements  indicate  that  the  vapor  is  very  nearly  fully  ionized.  The  high  velocities, 
in  particular,  are  too  large  ;o  be  accounted  for  by  any  reasonable  surface  tamper¬ 
ature  and  are,  in  themselves,  adequate  evidence  of  vapor  absorption.  Vapor  absorp¬ 
tion  ban  also  been  confirmed  by  direct  measurement  of  the  optical  opacity  of  the 
vapor . 

The  mass  density  experiment  was  intended  to  provide  additional  information  on  the 
details  of  the  interaction  process.  From  ti  qualitative  description  above,  one 
would  expect  that  the  mass  density  close  to  the  surface  would  show  a  maximum 
shortly  after  the  onset  of  the  laser  pulse,  tad  decrease  rapidly  -  the  rate  of 
evaporation  decreased.  The  decrease  in  the  rate  of  evaporation  v  ,d  be  expected 
if  the  direct  laser  radiation  was  prevented  from  reaching  the  surface  by  the  opaque 
'  apor.  A  mass  density  measurement  would  also  be  expected  to  show  the  presence  of 
shock  phenomena  in  the  vapor. 
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II.  EXPERIMENTAL  PROCEDURE 


A.  CONCEPT 

The  basic  concept  of  the  Mass  density  messurcL:->nt  is  to  direct  a  beam  of  high-energy 
electrons  (20-30  kev)  through  the  plasma  in  front  of  the  target,  and  to  measure  the 
attenuation  of  the  beam  as  a  function  of  time.  At  these  high  electron  energies, 
the  atoms  or  molecules  in  the  path  of  the  electron  beam  act  as  Rutherford  scatterers, 
with  the  bulk  of  the  scattering  being  caused  by  the  electrostatic  field  of  the 
nucleus.  Thus,  the  state  of  the  scattering  medium  is  not  of  great  significance, 
and  the  amount  of  scattt  mg  depends  primarily  on  the  number  of  atoms  in  the  path 
of  the  beam. 

Rutherford  scattering  is  highly  angular-dependant.  The  differential  scattering 
cross  section  for  an  atom  of  atomic  number  Z,  for  electrons  of  energy  £,  can  be 
written  as 


0(6)  = 


16  sin*  (6/2) 


On  insertion  of  numerical  values  and  electron  energies  in  kev,  this  equation 
becomes 


o(6)  =  1.29  X  10-^^  Z^/E~  sin^  (6/2)  cm^ 


The  total  cross  section,  determined  by  integrating  the  differential  cross  section, 
is  infinite.  In  actual  practice,  however,  one  will  collect  all  electrons  which 
are  scattered  by  less  than  some  angle  0^,  so  that  the  attenuation  cross  section  is 


=  1.29  X  (Z.'E)^  2n  I  sin  6  sin-*  (6/2)  d6 


f 

I  Sin 


a  =  1.29  X  (Z/E,  ^  4n  cot^  (6J2)  cm  ^ 


In  a  carefully  designed  experiment,  6^  can  be  made  as  small  as  10  radians,  so  that 


a  =  6.5x10-'^^  (Z/F.)^  cm^ 


For  graphite,  Z  ==  b,  and  a  typical  electron  beam  energy  is  25  kev.  Thus,  the 
attenuation  cross  section  per  particle  could  be  as  large  as  3.75  x  10“^^  cra^. 
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From  previous  work,  we  can  make  a  very  rough  estimate  of  the  mass  density  near  the 

target.  It  has  been  obser\'ed  that,  for  an  incident  energy  of  about  1  joule  on  10“^  ■ 

cm2  area,  about  10“^  grams  of  material  are  removed.  This  material  has  a  velocity 

in  the  vicinity  of  3  x  10^  cm/sec,  for  graphite.  If  we  assume  that  the  material 

leaves  the  target  over  a  time  approximately  that  of  the  laser  pulse  duration 

(5  X  10“®  seconds) ,  we  would  expect  a  mass  density  near  the  target  region  of 

about  6.7  X  10“^  gm/cm3. 

If  we  propagated  an  electron  beam  along  a  line  parallel  to  the  target,  through  the 
center  of  the  area,  the  beam  would  travel  a  distance  of  1.1  x  10~^  cm,  so  the  beam 
would  encounter  an  Integrated  density  of  about  7.4  x  10"^  gm/cm^,  or  3.7  x  10^® 
particles  for  graphite.  With  the  attenuation  cross  section  given  above,  the  heam 
would  be  attenuated  by  a  factor  of  e^®^,  so  that  it  is  obvious  that  a  much  lower 
coefficient  is  needed.  However,  lowering  of  the  coefficient  can  be  accomplished 
by  increasing  the  effective  aperture  of  the  detector  of  the  electron  beam.  By 
increasing  0^  to  5  x  10“2  radians,  the  attenuation  cross  section  is  decreased  by 
a  factor  of  25,  and  the  attenuation  is  now  a  factor  of  or  260.  This  is  still 

large,  but  it  is  in  the  vicinity  of  experimental  feasibility.  Other  effects,  such 
as  shielding  of  the  nucleus  by  outer  electrons,  would  be  expected  to  decrease  the 
attenuation  cross  section  to  some  extent.  Our  estimate  of  mass  density  may  also 
be  somewhat  high. 

The  experimental  observations,  then,  consist  of  determining  the  transmitted  electron 
beam  current  as  a  function  of  time.  The  time  is  to  be  related  to  the  time  of  onset 
of  the  laser  pulse.  The  integrated  density  encountered  by  the  beam  can  then  be 
determined  from  the  equation 


1  =  1^  exp  —  a 


where  I  is  the  instantaneous  beam  current,  the  current  in  the  absence  of  atten¬ 
uation,  a  a  calibration  constant  to  be  determined  either  experimentally  or  from  the 
equations  fo.  Rutherford  scattering,  and  p  the  mass  density. 

If  the  electron  beam  is  propagated  along  a  line  parallel  to  the  target  surface, 
at  a  fixed  distance  z  from  the  target,  and  several  observations  are  made  at  dif¬ 
ferent  locations  of  the  beam  from  the  central  axis  of  the  laser  beam,  it  is  pos¬ 
sible  to  determine  the  mass  density,  for  a  given  time,  as  a  function  of  i  and  the 
radius  from  the  central  axis,  provided  the  blowoff  pattern  is  s>’mmetrical  about 
this  axis.  The  mass  density  is  assumed  to  be  a  function  of  i,  r,  and  t.  The  ob¬ 
served  integrated  density  is  a  function  of  r,  y,  and  i.  The  integrated  density 
is  related  to  the  mass  density  by  the  equation 


c- 
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This  equation  is  one  form  of  Che  Abel  equation,  and  its  inverse  is 


A  number  of  numerical  methods  exist  for  the  determination  of  p  from  values  of  fp  dx 
determined  at  several  different  distances  y  from  the  axis. 

B.  APPARATUS 

Successful  completion  of  the  mass  density  experiment  presented  some  serious 
problems  earlier.  In  pa“ticular,  the  copious  amount  of  ultraviolet  radiation 
emitted  by  the  plasma  in  front  of  the  target  caused  some  serious  effects. 

This  radiation  caused  a  great  deal  of  photoelectric  emission  from  the  election- 
beam  detector  electrode,  which  in  the  initial  experiments,  was  a  Faraday  box. 

This  emission  caused  a  spurious  signal. 

The  detector  problem  was  eliminated  by  replacing  the  detector  with  a  scintillating 
material.  A  thin  metallic  film  on  the  surface  prevented  direct  electromagnetic 
radiation  from  reaching  the  scintillator.  A  photomultiplier  was  used  to  view  the 
light  emitted  when  electrons  passed  through  the  film  into  the  scintillator. 

A  second  problem  was  then  discovered;  the  ultraviolet  radiation  also  struck  the 
electrodes  of  the  electron  gun  which  produced  the  beam.  Since  these  electrodes 
were  at  a  high  negative  potential,  the  resultant  photoelectric  emission  caused  an 
appreciable  defocusing  of  the  beam,  and  again,  a  spurious  signal  was  observed. 

In  order  to  eliminate  this  problem,  as  part  of  the  current  work,  the  apparatus 
shown  schematically  in  Figure  1  was  designed.  Here  there  is  no  direct  optical 
path  from  the  target  to  the  electron  gun;  the  beam  passes  through  optical  baffles 
and  is  magnetically  deflected  to  pass  near  the  target. 

C.  CALIBRATION  OF  THE  ELECTRON  BEAM 

The  attenuation  coefficient  a  was  determined  by  passing  the  electron  beam  through 
a  thin  graphite  film.  This  coefficient  was  determined  to  be  about 

1.15  X  10~^  cm  ^/gm 


This  corresponds  to  an  attenuation  cross  section  of  2.3  x  10~^®  cm^,  which  is  the 
cross  section  to  be  expected  if  0^  =  0.04  radian. 

D.  EXPERIMENTAL  PROCEDURE 

Figure  2  shows  a  schematic  of  the  experimental  arrangement.  To  turn  on  the  electron 
beam  prior  to  the  actual  laser  pulse,  a  trigger  pulse  was  derived  from  the  fluores¬ 
cent  radiation  reflected  off  the  rotating  mirror  of  the  laser  oscillator.  The  posi¬ 
tion  of  the  pickup  of  this  radiation  could  be  adjusted  to  provide  a  trigger  about 
50  microseconds  before  the  laser  pulse. 
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LASER  BEAM 


Figure  1  ELECTRON  BEAM  APPARATUS 


The  laser  peak  power  and  pulse  shape  was  determined  with  the  aid  of  a  beam  splitter 
outside  the  experimental  chamber,  a  planar  diode,  and  a  Tektronix  519  oscilloscope. 
The  oscilloscope  was  triggered  Internally  with  the  laser  pulse  signal.  The  trigger 
pulse  from  the  oscilloscope  was  also  used  to  trigger  c.  second  oscilloscope  which 
recorded  the  actual  electron  beam  current.  A  second  beam  splitter  and  photodiode 
were  used  to  provide  a  tire-integrated  measurement  of  the  laser  output.  Both 
photodiodes  were  calibrated  against  a  water-filled  calorimeter  placed  in  the 
target  position  inside  the  apparatus. 

The  electron  beam  was  kept  in  essentially  the  same  position  during  all  of  the 
firings.  The  position  of  the  laser-beam  axis  with  respect  to  the  electron  beam 
was  adjusted  by  a  vertical  translation  of  the  lens  used  to  focus  the  beam  on  the 
target.  The  distance  between  target  and  beam  was  adjusted  by  a  translation  of  the 
target,  and  the  laser  beam  size,  at  the  target,  was  controlled  by  adjusting  the 
focusing  lens.  The  target  was  mounted  on  a  rod  which  could  be  rotated.  When  the 
laser  beam  was  deliberately  focused  off  the  axis  of  rotation,  it  became  possible 
to  rotate  the  target  between  firings,  presenting  a  fresh  surface  to  the  laser  beam, 
thus  preventing  excessive  hole-drilling  effects. 


III.  INTERPRETATION  OF  RESULTS 


A.  EXPERIMENTAL  RESULTS 

A  typical  e;rperitnental  result  is  shown  in  Figure  3,  where  the  integrated  density 
(the  mass  density  projected  along  the  line  of  the  electron  beam)  in  gm/cm^  is 
plotted  as  a  function  of  time  after  Initiation  of  the  laser  pulse. 

All  of  the  experimental  data  has  the  same  qualitative  appearance  —  a  peak  inte¬ 
grated  density  at  about  0.5  microsecond  after  the  laser  pulse,  and  observable 
effects  for  about  2  microseconds.  This  is  in  direct  contrast  to  the  results 
expected,  if  the  velocity  distributions  observed  with  the  microwave  cavity  tech¬ 
nique  dominated  the  density  close  to  the  target.  At  the  intensity  used,  the 
average  velocity  o.f  the  plasma  as  observed  with  the  cavity  techniaue  is  on  the 
order  of  10^  cm/sec.  Thus  the  peak  integrated  density  should  ocr  .  during  the 
laser  pui se  (in  less  than  3  x  10“^  seconds) . 

We  are  therefore  faced  with  the  serious  problem  of  accounting  for  the  observed 
electron  beam  results.  At  the  same  time,  it  must  be  kept  in  mind  that  the  ob¬ 
servations  with  the  microwave  cavity  indicate  much  higher  velocities  than  can  be 
inferred  from  the  electron  beam  results. 

The  microwave  cavity  data  indicate  that  the  average  velocity  of  the  plasma  at 
the  intensity  used  for  Figure  3  is  about  5  x  10^  cm/sec,  and  that  the  average 
kinetic  energy  is  about  4  x  10"^0  erg/atm.  The  total  mass  of  vapor  contributing 
to  the  observed  cavity  data  can  be  estimated  by  assuming  that  about  50  percent 
of  the  laser  energy  went  into  kinetic  energy  of  the  vapor.  For  graphite,  this 
results  in  a  mass  loss  of  about  0.25  x  10“®  gm/ joule  (for  the  example  given  in 
the  figure,  about  3.3  x  10“^  gm).  This  vapor  is  emitted  during  the  laser  pulse 
(about  50  nanoseconds),  so  that  the  mass  evolution  rate  is  about  6.5  gm/sec. 

For  the  sake  of  calculation,  let  us  assume  that  during  the  laser  pulse  the  vapor 
can  be  characterized  by  a  velocity  distribution  typical  of  evaporation  from  a 
surface  at  a  temperature  T.  From  kinetic  theory,  this  distribution  can  be 
written  as 

2 

F(v)dvd9.=  —  (il/2RT)~  cosO(exp-Mv^/2RT)dv  dU  (8) 
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where  F  (v)  dv  d9  is  the  fraction  emitted  in  the  velocity  range  v  to  v  +  dv ,  in 
the  direction  6,  and  the  solid  angle  rffi.  M  is  the  molecular  weight,  R  the  gas 
constant,  and  T  che  absolute  temperature.  From  this  distribution,  the  average 
kinetic  energy  is 

E  =  2  RT  erg/ mole  (9) 


•5 


and  the  average  velocity  is 
1/2 

3(2n  RT/M) 

V  =  -  cm/ sec  (10) 

4 

Fitting  the  kinetic  energy  with  the  observed  energy,  we  get  a  temperature  of 
1.45  X  10^°K,  and  from  the  velocity  we  get  T  =  lO^^K. 

Appendix  B  shovs  that  for  the  case  given,  the  integrated  density  can  be  expressed 
by  the  equation 

00 

1/2 

pdx  =  1.77  m  (M/2R  T)  '  z'-  \exp(-  Mz  ~/2R  Tt^)  -  exp  [-  Mz^/2R 

(11) 


where  m  is  the  mass  evaporation  rate  per  unit  area,  a  the  irradiated  radius,  z 
the  distance  of  the  beam  from  the  surface,  t'  the  laser  pulse  duration,  and  t  is 
the  time  after  the  start  of  the  laser  pulse.  In  plotting  the  expression,  the 
second  term  on  the  right-hand  side  is  zero  for  t  <  t'  .  We  can  evaluate  this 
expression  for  z  =  0.1  cm,  77  m  <?  =  6.5  gm/sec,  T  -  lO^^K,  )/  =  12  and  R  = 

8.31  x  10^  erg/mole-‘’K. 

p  (/*  =  9.85  X  10-^  I  exp  (-  7.22  x  )  -  exp  [-7.22  x  x  )-  ]  I 

(12) 

This  expression  is  also  plotted  in  Figure  3.  As  can  be  seen,  the  contribution 
for  the  high-velocity  species  is  significant  only  at  short  times,  and  we  must 
seek  another  means  of  accounting  for  the  data. 

This  discrepancy  between  the  microwave  cavity  data  and  the  electron  beam  data  is 
only  apparent  at  the  short  distances  of  the  electron  beam  experiment.  At  long 
distances,  the  electron  beam  experiment  is  not  applicable  (densities  would  be 
too  low  for  detection) ,  and  the  microwave  cavity  is  only  sensitive  to  the  elec¬ 
trons  present  in  the  high-velocity  plasma. 

In  order  to  account  for  the  observations,  it  appears  necessary  to  postulate  at 
least  two  components  of  the  vapor  emitted.  One  component  is  the  fast  vapor 
which  results  from  the  absorption  of  laser  energy  by  the  vapor  near  the  target. 
This  component  is  significant  only  at  "short"  times,  and  does  not  contribute  to 
the  bulk  of  the  ob ."arved  density.  The  second  component  can  be  assumed  to  be 
thermally  emitted  by  the  hot  target  surface,  both  during  and  after  the  laser 
pulse.  Most  of  the  vapor  emitted  during  the  laser  pulse,  however,  is  further 
heated,  and  is  therefore  not  a  major  concrlbutor  to  the  "slow"  component.  Let 
us  assume,  however,  that  perhaps  an  amount  of  vapor  equal  to  the  fast  component 
is  also  emitted  during  the  pulse,  but  is  not  accelerated,  and  therefore  has  a 
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temperature  characteristic  of  the  surface.  This  temperature  can  be  estimated 
from  the  mass  loss  rate,  by  means  of  the  approximate  equation 


m  =  7.2  X  10^  exp  (-  1.03  x  10^/T) 


(13) 


The  expression  is  an  approximation  to  the  published  data  on  the  rate  of  graphite 
evaporation  as  a  function  of  temperature,  and  is  more  fully  derived  in  Appendix  A. 
In  the  case  of  Figure  3, 


6.5 


and  tr  is  2  x  10  ^  cm^.  Hence  T  =  8.05  x  10^. 

Applying  this  temperature  to  the  equation  for  Integrated  density,  we  have 


pdx  =  1.0  X  10-‘^  t  exp(-  9.00  x  10-^'^/ )  -  exp  [-  9.00  x  lO'^^^/ft-S  x  10-^)^]  \ 

(14) 


This  expression  is  also  plotted  on  Figure  3.  With  allowance  for  a  longer  low- 
velocity  tail  on  the  high-velocity  component,  it  appears  that  the  combined  effects 
of  the  two  components,  emitted  during  the  laser  pulse,  can  account  for  the  early 
part  of  the  observed  integrated  density  function.  However,  the  temperature  which 
one  must  assume  in  order  to  fit  Equation  (11)  to  the  observed  integrated  density 
at  longer  times  is  only  about  lO^^K.  At  this  temperature,  the  rate  of  evapora¬ 
tion  is  negligible,  being  over  30  orders  of  magnitude  too  low.  We  therefore 
must  conclude  that  the  observations  can  not  be  reconciled  with  the  assumption 
that  all  of  the  evaporation  occurs  at  times  which  are  of  the  same  duration,  or 
perhaps  longer  by  a  factor  of  2,  as  the  laser  pulse.  We  must  postulate  that 
appreciable  evaporation  occurs  after  the  pulse. 

B.  SURFACE  COOLING 

First,  it  is  Important  to  determine  whether  or  not  such  evaporation  is  realistic. 
During  the  laser  pulse,  the  temperature  at  the  surface  is  nearly  constant,  since 
most  of  the  energy  incident  on  the  surface  from  the  laser  beam  is  used  to  eva¬ 
porate  material  from  the  surface.  If  all  of  the  incident  energy  is  absorbed  at 
the  surface,  and  if  the  thermal  conductivity  of  the  graphite  is  constant,  than 
the  surface  temperature  can  be  shown  to  be  approximately 


T  =  Tg  t  <  t' 

T  =  (2T^/n)sirr^  t7t  i  >  t'  (15) 
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The  decay  of  temperature  Is  much  too  fast  to  allow  appreciable  evaporation  from 
the  surface,  at  times  longer  than  the  laser  pulse.  There  are,  however,  two 
mechanisirs  which  may  allow  a  much  slower  decay  time.  In  the  first  place,  the 
thermal  conductivity  of  the  graphite  would  be  expected  to  be  much  higher  at  ele¬ 
vated  temperatures  than  at  room  temperature.  In  the  second  place,  there  is  a 
great  deal  of  thermal  reradiation  from  the  vapor  during  the  laser  pulse.  Since 
this  radiation  occurs  at  short  wave  lengths,  it  is  not  unreasonable  to  expect 
that  it  penetrates  relatively  deeply  into  the  surface,  and  that  consequently,  the 
surface  layer  is  heated  to  a  relatively  high  temperature  for  an  appreciable  depth. 
In  the  absence  of  losses,  the  surface  layer  will  remain  hot  for  an  appreciable 
time,  since  thermal  conductivity  into  the  interior  is  relatively  low  for  graphite. 
We  can  estimate  the  time  required  for  an  appreciable  change  in  surface  tempera¬ 
ture  by  assuming  that  the  temperature  of  the  graphite  was  at  the  end  of  the 
laser  pulse,  for  a  distance  S  into  the  surface.  If  only  thermal  conduction  is 
considered,  the  temperature  of  the  graphite  can  be  expressed  by  the  equation 

-n^  Di 

exp - - -  (16) 


where  D  is  the  thermal  diffusivity  of  the  graphite,  and  a  is  the  thickness  of 
the  graphite.  If  a  is  large,  the  temperature  at  the  surface  is 


T  =  ^2T/rrJ 


where  the  summation  has  been  replaced  by  an  integration.  The  integral  can  be 
shown  to  be  equivalent  to  the  integral 


7-  =  (nDt) 


u~^  (sinud)  (exp  -  DO  du 
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which  is  a  tabulated  function.  Figure  4  is  a  plot  of  T/T^  as  a  function  of 
2Dt/d^.  It  can  be  seen  that  the  decrease  in  temperature  is  negligible  for 
2Dt/S^  less  than  0.1.  Since  D  for  graphite  is  about  0.01,  the  time  required  for 
detectable  cooling  is  t  =  5.0  5^,  for  8  in  cm.  To  decrease  the  temperature  by 
10  percent  requires  20  8^  seconds. 

Thus,  if  the  heated  layer  were  only  3  microns  thick,  approximately  2  x  10  ^ 
seconds  would  be  required  to  cool  the  surface  by  10  percent,  by  conduction  alone. 
The  surface  temperature  is  not  expected  to  reach  more  chan  10,000®K.  A  3-mlcron 
layer  would  thus  store  approximately  6  joule/cm^  of  energy.  Since,  under  the 
conditions  of  the  experiment,  ti.tal  energy  deliveries  to  the  surface  are  on  the 
order  of  100  to  1000  joule/cm^,  it  appears  that  sufficient  energy  is  ava-*  '.able, 
even  if  only  a  fraction  penetrates  to  cause  huating  of  a  sufficiently  thick  layer. 

It  is  apparent,  then,  that  if  a  sufficiently  thick  surface  layer  is  heated,  the 
losses  by  evaporization  and  radiation  from  the  surface  must  dominate  the  cooling 
process.  The  loss  by  radiation  cannot  exceed  the  blackbody  limit 


I, 


'rad  =  ^  ^  wau/cm^. 


(20) 


The  loss  due  to  evaporation  can  be  approximated  by  the  expression 


'evap  =  ®  ^  10~'^  exp  (-1.03  x  10^/T )  watt/cm^ 


(21) 


Equation  (21)  is  also  derived  in  Appendix  B. 

Evaporation  losses  will  exceed  radiative  losses  for  all  temperatures  above  4200°K, 
and  therefore  radiation  can  be  neglected. 

For  times  during  which  the  effect  of  thermal  conduction  into  the  interior  can  be 
neglected,  the  cooling  rate  of  the  surface  can  be  approximated  by  the  cooling 
rate  that  the  surface  of  an  infinitely  thick  slab,  initially  at  a  temperature 
T^,  would  undergo. 

Suppose  we  have  a  very  thick  material  which  is  heated  w.  h  an  input  flux  I  at 
the  surface.  Then  the  temperature  at  a  time  t  ,  after  heating  for  an  incremental 
time  At',  is 


T 


I  Ht'/frr  p  CKt) 


1/2 


(22) 


If  /  is  a  function  of  t'  ,  then,  at  a  time  <  later  than  t'  ,  we  have 


dT/dt'  = 


(23) 


[n  pCK(t-t')] 


1/2 
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Since  heat  equation  solutions  can  be  superposed,  the  temperature  at  the  surface 
is  then  the  integral  of  this  expression  over  the  time  of  the  heat  flux  /  .  Al¬ 
ternatively,  if  /  is  a  function  of  temperature,  we  have 

.  - 1/2  -  1/2 

-[1(T)\-^  dT  =  (jipCK)  (i-f)  dt'  (24) 


which  is  to  be  integrated  from  t'  =  0  to  /'  •=  i.  Thus 


I  1/2  1/2 

KT)-^  dT  =  2  t  '  /  (npCK)  (25) 
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For  our  case,  I  i  T)  =  8x  10  exp  (-1.03  x  10  /  T  )  and  we  have 


It  1/2  1/2 

(exp  1.03  X  10^/T)dT  =  id  x  10^^  t  /(npCK) 


Equation  (26)  is  not  analytically  integrable.  However,  a  numerical  integration 
followed  by  a  curve-fit  results  in  a  fairly  good  approximation  to  the  integral 


3.81  X  10^  (exp  8.95  x  lO'^/T  -  exp  8.95  x 
=  1.6  X  iO^^  t^^^/ (n  pCK)^^~ 


(27) 


If  we  now  use  the  equation  for  m  as  a  function  of  T,  we  obtain  an  expression  for 
the  mass  evolution  rate  as  a  function  of  time: 


rf.  =  5.45  X  10-^  (t^''-+  1.89  X  10-^^  exp  8.95  x  lO'^/T^r^-^^ 


(28) 


Over  most  of  the  range  in  t  ,  the  evolution  is  reasonably  well  represented  by  the 
simpler  expression 


=  5.45  X  10  ^-0.575  gm/cm^  -sec 


(28a) 
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Equation  (27)  can  also  be  used  to  provide  an  approximate  expression  for  the  tem¬ 
perature. 


i  <} 


T  =  8.95  X  I0‘*  (In  5.29  x  10^^  t^^-) 


(27a) 


Equation  (27a)  is  plotted  in  Figure  5.  As  can  be  seen,  the  temperature  is  very 
nearly  constant  over  a  wide  range  in  time. 

C.  ESTIMATED  DENSITY 


We  have  shown  that  it  is  not  unreasonable  to  expect  that  the  temperature  of  the 
graphite  surface  will  remain  high  fo’.  a  relatively  long  time  following  the  laser 
pulse.  To  be  sure,  the  temperature  drops  rapidly  at  first,  but  it  remains  in 
the  vicinity  of  5000‘’K  for  an  appreciable  time. 


We  will  use  the  same  form  used  to  derive  the  density  for  a  constant  temperature 
evolution.  The  integrated  density,  due  to  evaporation  during  a  short  increment 
of  time  dt'  ,  can  be  expressed  as 


M 

2  g  2RT(t-t')^ 


dt' 


(29) 


Although  both  m  and  T  in  this  expression  are  functions  of  time,  the  value  of  T 
varies  very  slowly  in  time.  Thus  we  are  justified  in  replacing  T  by  an  average. 
The  average  value  of  T  can  be  estimated  by  using  the  value  of  T  that  would  give 
the  same  total  mass  emission  in  a  time  of  about  1  microsecond.  The  total  mass 
emission  in  10“^  seconds,  according  to  the  approximate  expression  for  mass 
evaporation  rate,  is  3.6  x  10“^  gm/cm^.  This  would  correspond  to  a  constant 
temperature  of  5.25  x  lO^'K.  The  integrated  density,  at  a  time  t  ,  then  becomes 


J!j 


p  dx 


2  X  5.45  X  10-“^  A 


V  ” 


e  2RT(t-t’J 
(i  -  t')^  (i')0.575 


dt' 


(30) 


If  we  let 


9 


u 


2  HT  r/  M  z*‘  ,  we  have 


(2)  (5.45  X  lO-'h  A  /  M  \  0.2i25 
Z  \  2«7'/ 

\:  TT  '  ' 


1 


e 

(u')0.575  (^_^')3 


du' 

(31) 
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Figure  5  COOLING  OF  GRAPHITE  BY  EVAPORATION 


where 

u  =  2J  y.  10^  t 


The  integral  again  is  not  analytic,  but  can  be  numerically  evaluated.  The  inte¬ 
gral  is  given  in  Figure  6. 


The  maximuia  value  of  the  integral  is  0.68,  occurring  at  u  =  1.2.  For  the  case 
in  Figure  3,  A  =  2  x  10“^  cm^,  giving  a  maximum  integrated  density  of  1.5^  x  10"^ 
gm/cm^  at  a  time  of  4.4  x  10“'  seconds.  The  time  is  in  good  agreement  with  ob¬ 
servation,  but  the  value  of  the  maximum  is  about  two  orders  of  magnitude  too  low. 
However,  it  should  be  recalled  that  the  room  temperature  value  of  graphite  ther¬ 
mal  conductivity  was  used  to  calculate  the  mass  loss  rate.  If  the  thermal  con¬ 
ductivity  were  to  be  increased  by  a  factor  F,  the  equation  for  mass  loss  rate 
would  be  increased  by  the  factor  F 0-575^  Furthermore,  it  is  not  at  all  unreason¬ 
able  to  postulate  that  the  area  heated  by  re-radiation  from  the  vapor  would  be 
larger  than  the  focal  area  of  the  beam.  An  area  increase  by  a  factor  of  10  does 
not  appear  to  be  too  much  for  a  focal  diameter  of  only  0.05  cm.  The  high-velocity 
plasma  travels  a  distance  of  about  0.25  cm  during  the  laser  pulse  so  that  irradi¬ 
ation  of  a  comparable  diameter  at  a  high  intensity  level  might  be  expected  to 
occur.  An  increase  in  thermal  conductivity  by  a  large  factor  might  alsc  be 
realistic.  It  is  therefore  not  unreasonable  to  normalize  the  calculated  inte¬ 
grated  density  to  fit  ti—  observations.  The  normalization  is  shown  in  Figure  7, 
along  with  the  experimental  observations. 


One  of  the  disturbing  factors  in  the  preceding  development  is  that  it  was  neces¬ 
sary  to  postulate  chat  the  temperature  and  area  vere  higher  than  those  initially 
assumed.  Also,  the  calculated  Integrated  density  does  not  fall  off  as  rapidly, 
at  long  times,  as  that  observed.  This  latter  is  of  course  a  result  of  the 
neglect  of  the  effect  of  thermal  conduction  into  the  interior. 


Let  us  consider  an  alternate  approach  to  the  problem.  Suppose  we  assume  that  a 
finite  thickness  of  material  is  heated  to  the  initial  temperature  and  that  the 
thermal  conductivity  of  this  very  hot  layer  is  sufficient  so  that  a  constant 
temperature  is  maintained  throughout  the  layer,  then  the  rate  of  cooling  (ne¬ 
glecting  thermal  conductivity  into  the  interior)  is  given  by 


Figure  6  INTEGRATION  OF  EQUATION  (32) 


dT 

IT 


8>.I0l3(exp-^-03^l0^/r) 

pCS 


exp  1.03  X  I(^/T  dT  = 


8  X  10^^  t 
pC8 


(33) 


The  expression  on  the  left  has  already  been  evaluated,  and  the  resulting  expres¬ 
sions  for  mass  evaporation  rate  and  temperature  are,  using  pC  =  5  joule/cm^, 


1.02  X  10-^ 

(,/t  *  2.!«  X  io-'%  «•«  ■<  /•'" 

T  =  8.95  X  10'^  [In  (4.9  -  10^  t/S  +  exp  8.9  x 


(34) 


(35) 


Equation  (35)  indicates  a  slow  variation  in  T  for  most  of  the  elapsed  time. 
Again,  we  can  estimate  an  average  temperature  by  choosing  a  constant  temperature 
which  gives  the  same  emission.  In  one  microsecond.  Equation  (34)  gives  an 
emission  of  1.0  x  10"^  gm/cm2,  for  an  average  temperature  of  6320°K,  for  a  S  of 
3  X  10“^  cm,  and  T  -  10^®K.  We  must  then  evaluate  the  equation 


f 


p  dx  ^  -  (1.02  X  W 

VT 


t/»- 


2RT(t-t’r 


9/  T*  1.13 


(‘-’’r  (1/8  +  2.30  X  /O-^"  X  ^OVT'o; 


T  d\' 


(36) 

Again,  this  equation  must  be  numerically  integrated.  If  we  let 


2RT  t 


V2RT 

“IT 


we  have 


/ 


p  dx 


2  A 


x/ir 


(1.02  X  10 


/  ,,  \  -0.075  /•“ 


(u-u')^  (b+u')‘-^^ 
(37) 


du ' 
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u  =  2.97  X  10^  t, 
b  =  1.6  X  10-^ , 

and 

I  p  dx  =  9.52  X  10-‘^  A  I 

Figure  8  gives  the  results  of  the  numerical  evaluation.  The  integral  has  a 
maximum  value  of  3.85  occurring  at  »  =  0.92.  This  corresponds  to  an  integrated 
density  (using  A  -  2  x  10~^)  of  7.3  x  10"^  at  a  time  of  3.1  x  10"7  seconds.  The 
integrated  density  is  about  a  factor  of  6  too  low,  however  it  is  again  reasonable 
to  assume  that  the  area  evaporating  is  larger  than  the  Irradiated  area. 

Figure  7  shows  the  results  of  these  two  approaches,  as  well  as  the  observed  re¬ 
sult.  Both  calculated  curves  are  normalized  to  give  the  same  peak  integrated 
density. 

The  curve  for  the  infinite  slab  case  appears  to  provide  a  better  approximation 
at  times  near  the  peak  of  the  integrated  density,  whereas  the  curve  for  the 
finite  slab  appears  best  for  much  longer  times.  At  short  times,  the  two  curves 
are  approximately  equal  (after  normalization),  and  neither  accounts  for  the  ob¬ 
served  density.  The  dashed  line  shows  the  difference  between  the  calculated 
and  observed  values  at  short  times.  This  difference  is  assumed  to  be  the  high- 
velocity  component  of  the  vapor  emitted  by  the  target,  i.e.,  the  vapor  which 
has  been  accelerated  by  gaining  energy  from  the  laser  beam. 

D.  OBSERVATIONS  AT  LARGER  DISTANCES 

The  data  presented  in  Figure  3  (and  repeated  in  Figure  7)  are  typical  of  all  the 
data  taken  at  short  distances  from  the  target.  To  be  sure,  there  are  variations, 
but  not  of  a  significant  nature.  Figure  9  presents  some  data  obtained  at  a 
longer  distance  (3mm  from  target).  Again,  this  doer  •»ot  differ  greatly  from 
that  obtained  at  the  shorter  distance.  There  is,  however,  a  tendency  towards  a 
separation  of  the  high-  and  Itw-velocity  components.  This  tendency  is  presumably 
fortuitous,  since  It  is  not  always  observed.  It  can  be  accounted  for,  if  the 
relative  amount  of  the  higher  velocity  component  is  increased  over  the  one  given 
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Figure  9  INfEGRATED  DENSITY  VERSUS  TIME 
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IV.  CONCLUSION 


The  use  of  the  electron  beam  technique  to  determine  the  mass  density  of  the  laser 
blowoff  has  turned  out  to  be  of  doubtful  value,  not  because  the  experiment  was 
unsuccessful,  but  because  the  technique  responded  to  effects  which  were  not 
properly  a  part  of  the  high  Intensity  interaction. 

It  was  originally  hoped  chat  Che  mass  density  approach  would  confirm  the  details 
of  theoretical  models  on  the  interaction.  However,  the  experimental  results 
themselves  indicate  that  the  bulk  of  the  experimental  results  are  due,  not  to 
nlasma  created  during  the  time  duration  of  the  laser  pulse,  but  to  the  vapor 
which  has  been  thermally  emitted  from  the  hot  target  after  the  pulse. 

This  result  is  in  itself  of  interest,  since  it  can  account  for  phenomena  which 
’  ive  been  observed  in  other  experiments.  One  such  phenomena  is  the  observation 
of  radiative  decay  of  species  driven  from  a  solid  target  by  laser  radiation.* 

The  observations  tend  Co  give  lifetimes  longer  than  could  reasonably  be  obtained 
if  the  vapor  left  the  vicinity  of  the  target  in  a  time  which  was  characteristic 
of  thermal  evaporation  during  the  laser  pulse  alone. 

In  order  to  use  Che  electron  beam  technique  to  probe  events  that  occur  close  to 
the  surface,  it  is  necessary  to  concentrate  on  the  first  event  during  the  laser 
pulse.  In  the  experiments  that  have  been  performed,  the  whole  event  was  observed, 
so  that  the  oscilloscope  traces  of  the  electron  beam  current  were  observed  for 
times  of  typically  2  to  4  microseconds;  consequently,  when  the  effort  to  observe 
the  whole  event  was  made,  the  time  resolution  necessary  to  perform  an  analyses 
0  for  the  first  few  tens  of  nanoseconds  was  lost. 

We  have  now  found,  however.  Chat  the  whole  event  is  not  meaningful  for  the  high- 
intensity  interaction.  In  fact,  the  only  part  of  the  high-intensity  interaction 
which  can  be  observed  with  the  electron  beam  technique  is  the  part  that  occurs 
before  the  leading  edge  of  the  low-velocity  vapor,  thermally  emitted  from  the 
target,  has  reached  the  position  of  the  electron  beam. 


•  Private  communication  from  T.  Wentink,  Avco/SSD. 
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APPENDIX  A 


EVAPORATIVE  PROPERTIES  OF  GRAPHITE 

The  vapor  pressure  of  graphite  is  published  as  a  function  of  temperature  e.*  If 
the  vapor  pressure  is  known,  the  mass  evaporation  rate  can  be  determined  from 
kinetic  theory  to  be 

m  =  p/(2:tRT)^^^  (A1) 

where  p  is  the  pressure  in  dyne/cm^,  and  R  is  the  gas  constant  per  gram. 

From  the  published  data,  we  then  obtain: 

p  =  1  torr  T  =  3859®  K  m  =  3.26  x  10-3 
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This  can  be  fitted  quite  accurately  with  the  equation 

m  =  1,2  X  10^  exp  (— 1.03  X  lO^/T) gm/cm^-sec  (A2) 

The  latent  heat  of  evaporation  of  graphite  is  172  kcal/mole,  6.0  x  10^  joule/gm. 
Since  the  energy  required  to  evaporate  1  gm  of  graphite  at  a  temperature  T  is 
£'  =  6  X  10^  +  2RT/M,  the  power  per  unit  area  involved  in  evaporation  is 

I  =  (6  X  I0‘^  +  1,39  T)  m  watt/cm ~  (A3) 

This  can  be  approximated  by 

/  =  8.0x10^^  exp  (-1.03x10^/7)  wall/cm^ 


Handbook  of  Chemisiry  and  Physics,  41st  Edition. 


APPENDIX  B 


DERIVATION  OF  DENSITY  EQUATIONS 


Let  us  estimate  the  density  to  be  expected  from  the  vapor  evaporated  by  a  hot 
surface.  To  perform  the  calculations,  we  will  assume  that  the  vapor  is  Maxwellian, 
with  a  temperature  T  characteristic  of  the  surface.  From  kinetic  theory,  the 
velocity  distributio'i  function  of  a  gas  evaporating  from  a  surface  at  a  temperature 
T  is 

,V (v)  dvdn  =  (2/77)  (M/2R T)-  cos  e  (exp  -  .»/ v^/2 RT)  dv  dCl  (B1 ) 


and  is  the  fraction  emitted  in  the  direction  6  ,  measured  from  the  normal,  in  the 
solid  angle  dQ  and  in  the  velocity  range  from  v  to  v  +  dv  .  If  the  total  mass 
evaporation  rate  per  unit  area  is  m,  then  the  mass  evaporation  rate  per  unit  sur¬ 
face  area,  in  the  same  velocity  range,  direction,  and  solid  angle,  is 

m  N(v)  dv  dQ 


We  can  determine  the  density  at  a  distance  r  from  the  point  of  evaporation,  due 
to  the  mass  emitted  in  a  time  dt'  ,  by  first  determining  the  mass  evaporated  and 
dividing  this  by  the  volume  occupied  by  the  vapor.  The  evaporated  mass  in  the 
velocity  range  v  to  v  +  dv  ,  direction  0,  solid  angle  dfi,  and  timed  (/«',  per  unit 
surface  area  is, 

dm  =  mi\(v)dvdQ  dt'  (B2) 

The  volume  occupied  by  this  vapor  is  a  shell  of  area  r~ dil  and  thickness  v  dt' 
Thus  the  density  is 

dp  =  m(2/n)(M/2RT)  [v^  cos  e/r^  exp(-Mv^/2RT)]dv  (B3) 

At  a  time  t  ,  the  velocity  v  corresponding  to  evaporation  at  a  time  t  s 


V  =  r/(t  —  t')  (B4) 


where  t'  is  the  time  of  evaporation  r  i  t  is  the  time  of  observation. 
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Thus 


dv  =  rdt'/(l-t'J^ 


(B5) 


Also,  cos  0  =  z/r  ,  where  i  is  the  distance  normal  to  the  surface.  We  then  have 


dp  =  m(2/7T)(.\l/2RT)-  z  \cxp[-.\ir-/2RT(t-t')-]\dr 


(B6) 


The  experimental  observations,  however,  are  of  the  integrated  density  along  a  line 
parallel  with  the  surface.  Equation  (B6)  can  be  integrated  along  a  line  parallel 
to  the  surface.  Letting  r-  =  x~  +  +  z~  ,  have 
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m  z  (2/7T )  (\f/2RT)- 


exp 


~-M(yKz^)~ 

2RT(t-t')~ 


dx 
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m  z(2/\fn) 


exp 


2RT(i-t')^ 
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(B7) 


This  expression  can  be  integrated  over  t'  ,  if  desired.  However,  we  have  not 
considered  the  effect  of  the  finite  source  area.  If  we  consider  a  circular 
source  area,  we  have,  for  the  density  integrated  along  a  line  passing  through 
the  normal  to  the  center  of  the  source. 
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p  d\ 


M »  - 


JRT(t-i’)- 


dt’ 


The  integral  is  equivalent  to  the  integral 


(B8) 


(B9) 


-31- 


where 


F  gure  Bi  INTEGRATION  OF  EQUATION  (B9) 


Even  i£  z  is  less  than  a,  it  can  be  seen  from  Figure  B1  that  the  integral  over 
the  surface  area  is  not  a  very  rapidly  varying  function.  It  therefore  appears 
reasonable  to  use  the  value  of  the  integral  of  the  point  at  which 


is  a  maximum  The  maximum  of  Equation  (B8)  occurs  approximately  when 


1.5  = 


or  when 


2RT(t-t')‘ 


,  c  2  2 
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When  h  is  less  than  5,  the  graph  of  Figure  B1  is  applicable, 
the  integral  can  be  approximated  quite  well  by 


l^hen  b  is  larger, 


pdx  =  3.28  maz  ' 


(t  —  t')  ^  exp 


2RT(t-t')^J 


when  a'>  2z  ,  and 


d  I  p  dx  =  3.o4  m  a- z  [-—]  (t-t  )  ^  exp  -  - - 

J  \2«//  (_  2RT(t-t'r_ 


when  a  <  z 
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